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Abstract: 

Using a semi-classical treatment of the linear a model, we simulate the dynamical evolution of an 
initially hot cylindrical rod endowed with a longitudinal Bjorken scaling expansion (a "Bjorken 
rod"). The field equation is propagated until full decoupling has occurred and the asymptotic 
many-body state of free pions is then obtained by a suitable Fourier decomposition of the field 
and a subsequent stochastic determination of the number of quanta in each elementary mode. The 
resulting transverse pion spectrum exhibits visible enhancements below 200 MeV due to the para- 
metric amplification caused by the oscillatory relaxation of the chiral order parameter. Ensembles 
of such final states are subjected to various event-by-event analyses. The factorial moments of the 
multiplicity distribution suggest that the soft pions are non-statistical. Furthermore, their emis- 
sion patterns exhibit azimuthal correlations that have a bearing on the domain size in the source. 
Finally, the distribution of the neutral pion fraction shows a significant broadening for the soft 
pions which grows steadily as the number of azimuthal segments is increased. All of these features 
are indicative of disoriented chiral condensates and it may be interesting to apply similar analyses 
to actual data from high-energy nuclear collision experiments. 
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I. INTRODUCTION 



The possibility of forming disoriented chiral condensates (DCC) in high-energy nuclear 
collisions has gained significant attention in recent years as a means for testing our under- 
standing of chiral symmetry restoration A crossover from the normal phase, in which 
chiral symmetry is broken, to a phase with approximately restored chiral symmetry is ex- 
pected to occur at temperatures of a few hundred MeV. Such energy densities are readily 
reached during the early stage of an ultrarelativistic nuclear collision and one may therefore 
expect that approximate restoration of chiral symmetry occurs transiently in the hot col- 
lision zone. The subsequent relaxation towards the normal vacuum has a non-equilibrium 
character due to the very rapid expansion of the system. As a result, the pion field acquires 
large-amplitude long-wavelength oscillations which may have observable consequences. The 
suggested signals include an excess of isospin-aligned soft pions with an associated anoma- 
lously broad distribution of the neutral pion fraction and, in the electromagnetic sector, 
a significant excess of photons |IJ and dileptons |11| , [12| . 



A few experimental DCC searches have already been carried out [0,0] but no signals 
were discerned so far. On the theoretical side, a considerable number of explorative studies 
of this hypothetical phenomenon have been made over the past several years. However, 
due to the complicated nature of the problem, most studies have treated fairly idealized 
scenarios. Therfore the practical observability of the phenomenon is not yet clarified. The 
purpose of the present work is to treat a somewhat more refined scenario and thus obtain 
a better basis for making an assessment of the prospects for observing DCC signals and, in 
the process, provide some guidance for the analysis of the experimental data. 

The most popular tool for DCC studies has been the SU(2) linear a model which describes 
the 0(4) chiral field <p = (cr, 7r) by means of a simple effective quartic interaction. This 
framework has provided instructive insight into the features of the non-equilibrium DCC 



dynamics [15-33H- The present study is based on a semi-classical treatment of the standard 
SU(2) linear a model [E3]. The dynamical variable is thus the 0(4) chiral field <p(r) = (<r, n) 
which is subject to a self-interaction of the form V = \[4> 2 — v 2 ] 2 — Ha and, consequently, 
it is governed by a second-order non-linear equation of motion, 

[□ + A(0 2 -^ 2 )]0 = He a , (1) 

where e a denotes a unit vector in the a direction. (Here and throughout we often omit 
factors of % and c in order not to clutter the formulas.) 

In the present study, we wish to consider idealized scenarios that exhibits some of the 
most important features expected in real collision events, namely rapid longitudinal expan- 
sion and finite transverse extension. Generally, the numerical propagation of fields exhibiting 
significant flow patterns, such as rapid expansion, is practically difficult due to the phase 
oscillations caused by the local boost. However, for idealized scaling expansions this compli- 
cation can be eliminated by suitable variable transformations. For the longitudinal scaling 
expansion considered here pi] , it is convenient to replace the usual fixed-frame space-time 



variables (x,y,z,t) with the comoving variables (x,y,rj,r), 

t = t cosh f] , z = rsinhr/ . (2) 

Thus, t = (t 2 - z 2 ) l/2 is the proper time experienced in a system boosted along the z axis 
with the local rapidity y equal to the value of n = |ln(|^f). The corresponding form of 
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the field equation of motion is then obtained from the transformation of the d'Alembert 
operator, 



□ ee dl-dl-dl-dl = \d T rd T -dl-dl-^dl. (3) 

The equation of motion can be readily solved numerically by application of the leapfrog 
method, once the initial field 4>(r) and its time derivative ij)(r) are specified. (To help avoid 
confusion between the transverse coordinate y and the rapidity y, we shall preferentially 
denote the position in the transverse plane by p = (x, y) in the following, in analogy with 
the use of k = (k x , k y ) for the transverse wave number.) 



II. BJORKEN MATTER 

We wish to first consider a macroscopically uniform system, corresponding to infinite 
matter undergoing a longitudinal Bjorken scaling expansion. Without the scaling expansion, 
such a system would model infinite uniform matter. We refer to the present expanding 
generalization as Bjorken matter to capture the essential character of the system by a simple 
term. Just as ordinary matter approximates conditions prevailing in the bulk region of a 
finite system, such as a large atomic nucleus, Bjorken matter approximates the conditions 
in the interior of the rapidly stretching firestreak generated in a central high-energy collision 
of two large nuclei. 

For the numerical calculations, we represent the field on a cartesian lattice in space and 
impose periodic boundary conditions in all three directions. The initial field configuration 
is prepared by means of the sampling method developed in Ref. [p5|1 . Thus an initial tem- 
perature To is specified and a field 4>'(r) and its time derivative tp'(r) are sampled from 
the corresponding thermal ensemble. This field configuration represents a non-expanding 
system in thermal equilibrium. In order to obtain an initial field configuration suitable for 
the expansion scenario, the coordinate z is subsequently identified with the longitudinal 
variable rj, 

(f)(p,T],To) = (f>'(p,Z = 7]) , (4) 

i>(p,V,r ) ee if)'{p,z = rj) , (5) 

where we use r = 1 fm/c for the initial value of the proper time, in accordance with common 
practice. In this manner it is assured that the local environment at (p,rj), when analyzed 
in a frame boosted with a rapidity y equal to rj, approximates that of matter in equilibrium 
at the specified temperature T . 



A. Dynamical evolution 

It is instructive to note the evolution of some of the key quantities. Most important is 
the time dependence of the average field strength, <fi = ((f)), and the associated dispersion 
of the field fluctuations around that average, 5(p = [(<fi 2 ) — (po] 1 ^ 2 , where (fio ee |0|. The 
former represents the 0(4) chiral order parameter, while the latter reflects the presence of 
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quasiparticle excitations relative to that constant background field. The dynamical evolution 
of these quantities is displayed in Fig. [l] for a range of initial temperatures. Initially, at 
t = To, the system is hot and the order parameter is significantly reduced relative to its 
vacuum value of f n = 92 MeV, with the field fluctuations being substantial (and the more 
so the higher the value of To). The longitudinal expansion effectively cools the system, 
causing the field fluctuations to subside (the field dispersion 5(p evolves approximately as 
1/\/t at large times). This in turn makes the effective potential gradually revert to its 
vacuum form and the order parameter exhibits a corresponding relaxation. The specific 
oscillatory evolution is a result of the delicate balance between the expansion rate and 
oscillation frequency of the evolving effective potential. In the intermediate temperature 
range where the phase cross over occurs the order parameter exhibits a significant overshoot. 
These oscillations in the order parameter are reflected in the effective mass tensor for the 
quasiparticle excitations [^] and, accordingly, the field fluctuations (which depend inversely 
on the effective mass) exhibit corresponding undulations superimposed on the overall steady 
decrease in time. 

It is helpful to present the combined information in Fig. |as a projection onto the chiral 
phase diagram, in which the order parameter is the abscissa and the field fluctuation is the 
ordinate. The former characterizes the medium-modified vacuum and the latter is a measure 
of the degree of excitation relative to that vacuum. The result is displayed in Fig. |2|. This 
figure includes the curve connecting the equilibrium points, starting from <p gs = (f n , 0) 
at T = and approaching the vertical axis at large temperatures as chiral symmetry is 
gradually approached. This curve represents the adiabatic path along which the system 
would evolve if the expansion were sufficiently slow. We note that the scaling expansion 
quickly brings the system away from equilibrium, as the cooling of the field fluctuations is 
immediate while the order parameter requires some time to adjust to the changed effective 
potential. This inertia leaves the order parameter on the inside of the minimum in the 
effective potential and it therefore overshoots towards the outside once it gets moving. This 
pattern is repeated subsequently but on a smaller scale. As a result, the evolution is fairly 
soon confined to within a narrow neighborhood around the adiabatic path. Nevertheless, the 
oscillations in the order parameter persist for a long time with a slowly decreasing amplitude. 
Included on the display is also the critical boundary within which the field is supercritical 
and spontaneous pair creation occurs |25j. We note that the dynamical paths stay safely 



away from this boundary, regardless of the value of To, so that no signals of supercriticality 
should be expected for purely longitudinal expansions, as has already been noted in earlier 
studies 

In order to get an idea of what outcome might be expected, it is helpful to consider the 
time evolution of the effective pion mass, 

/4(t) = X[4 > 2 + 64 > 2 + 26n 2 -v 2 ] , (6) 

where 5(p 2 = 5a 2 + 5ti 2 + + 5n 2 is the total field fluctuation and 5n 2 denotes the contribu- 
tion from the field component in the particular isospin direction of the pion considered. This 
quantity is easily extracted as a by-product of the dynamical calculation and it is displayed 
in Fig. |3[ Since the phase trajectory stays away from the critical region (see Fig. |2|), fi 2 re- 
mains positive throughout. Its evolution can be characterized as an overall steady decrease 
(approximately inversely proportional to r) superimposed by a fairly regular oscillation. 
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It is often instructive to discuss the physics within the framework of the mean-field 
approximation, which describes the system as a gas of independent quasi-particles endowed 
with the time-dependent effective mass obtained self-consistently as indicated above (^|). In 
particular, one may consider the evolution of a single pion mode k being subjected to the 
calculated effective mass ^(t). This problem was studied recently in Ref. |35|] and it was 
found that the result can be conveniently expressed by means of the number enhancement 
coefficient X^. This quantity expresses by how much the occupancy of the given mode is 
enhanced as a consequence of the prescribed time dependence of its mass. Figure |4] shows 
the resulting enhancement coefficient as a function of the pion momentum, obtained for the 
four cases depicted in Fig. 

On the basis of this schematic analysis, one would expect a significant enhancement for 
the soft pion modes. It appears that the effect is largest for the intermediate values of T , 
as might have already been expected from the behavior seen in Fig. |[ Moreover, a careful 
inspection of the results reveals that Xn(p) has two components: a steady decrease resulting 
from the overall decay of fi n (t) and a bump caused by the parametric amplification of modes 
matching the frequency of the oscillatory part of fi n (t), with the former being by far the 
dominant component. (The location at the bump at 200 — 250 MeV/c corresponds to a pion 
energy of about half the a mass, as should be expected for the parametric resonance effect 



B. Analysis of the final state 

As noted above, the continual longitudinal expansion causes the local field fluctuations 
to subside steadily in the course of time. For the idealized one-dimensional expansion char- 
acteristic of Bjorken matter, 5<p 2 decreases as 1/r for large times (and faster than that when 
a system with a finite transverse cross section is considered, as in the next section, due to the 
additional dilution caused by the transverse expansion). Therefore, after a sufficiently long 
time the system will have decoupled into independently evolving modes of ever decreasing 
field amplitude. In this asymptotic regime, it is possible to extract well defined values of the 
observables by suitable analysis of the field configuration. It is useful to note that for large 
times the longitudinal velocity of a given part of the system is given by v z — > z/t = tanh(r/). 
Thus, in that limit, the coordinate rj equals the rapidity y = tanh _1 (t> 2 ). 

When analyzing the asymptotic field in the Bjorken expansion scenario, it is natural to 
perform a Fourier transformation in the transverse plane, 

tfkfo) = / J£ <KP,r,) e** , (7) 

= J |f MM) e " Jk p . (8) 

where k = (k x , k y ) denotes the transverse wave vector and Q± = L x L y denotes the transverse 
area of the spatial lattice. In the present study, we are only interested in observables 
based on the three pion components of the chiral field. Moreover, while the dynamical 
treatment of these components is carried out in a cartesian representation, ir = (tti, tt 2 , tt-s), 
the observations are made in the spherical representation, 7r = (tv_, 7t , 7r + ), where the 
individual components represent the charge states of the pion, 
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1 . . 

7T ± = -j= [7Ti ± m 2 \ 



7T 



V"3 • 



(9) 



The present treatment yields the time evolution of the classical pion field, 7r(x, y, r), r), as 
described above. In order to make contact with the physical observation, it is necessary to 
transcribe the emerging field into a quantal many-body state. For this purpose, we associate 
the asymptotic free field with a coherent state given by 



x) = exp \ / dr l \x\M) ■ at(rj) - xUv) ■ a M] \ |0) 



(10) 
(11) 



Here a^'(p, rj) annihilates a pion of type j located at (p, rj). It satisfies the usual commuta- 
tion relation, 



exp I d 2 p dr] x(p, v) • at (P> v) ~ X*{P, v) • «(P> V) \ 1°) • 



[a^(p',r, r ),a^(p, V y\ 



5 fj 8^(p'-p) 6{rf- V ). 



(12) 



Correspondingly, (77) annihilates a pion of type j located longitudinally at rj and having 

the transverse wave vector k, with [a^, (rf), (vV] = ^'j^u.'^iv' 

The coefficients in the expression (|10D for the coherent state are given simply in terms 
of the Fourier expansion coefficients in Eqs. (0-§). 



xM 



2 k y/2m k 



(13) 



where the transverse mass is given by m\ = m 2 + k 2 + ky. The corresponding Fourier 
expansion yields the complex isovector field appearing in the representation (|TT1) of the 
coherent state, 



x(p,v) 



ExM 



(14) 



-L k 



It presents a convenient encoding of the pion field configuration, i, e. the local field strength 
n(p,rj) and its time derivative. We note that the coefficient Xk(v) * s represented by the 
basic annihilation operator a^(rj), while %(p, rj) is represented by a(p, r/). 

The coherent state \x) has a number of convenient properties. Most importantly, it is 
an eigenstate of the annihilation operator, 



aM\x) = Xuiv)\x) , 
a(p,r])\x) = x(p,v)\x) 



(15) 
(16) 



with the eigenvalues being the corresponding x coefficients. The two-point density matrix 
for the state \x) is then given simply in terms of the complex x field, 



(x\a\p,v)a(p',r]')\x) = X*(p,v) x(p',v') ■ 



(17) 



Thus the local density of pions with isospin component j is given by the corresponding 
diagonal element, d?N^ / d 2 pdrj = |x^(p,?7)| 2 . Furthermore, the mean number of such 
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pions emerging with transverse wave vector k and a rapidity y in the interval (yi, y?) is then 
given by 

4 J) (yi,y 2 ) = [ y2 dv\x^(v)\ 2 , (is) 

where we have utilized the asymptotic equivalence of the longitudinal coordinate rj and the 
rapidity y. 

In the numerical treatment, the cartesian lattice has spacings of 0.2 fm in the transverse 
plane and Arj = 0.2 longitudinally. Thus the numerical resolution in the longitudinal direc- 
tion corresponds to a rapidity slice of width Ay = 0.2. In the analysis of the results, we 
consider rapidity intervals (yi, y2) that have unit length, corresponding to a physical rapidity 
resolution of one. The corresponding sources are thus made up of five contiguous rapidity 
slices and we shall refer to them as lumps. 



C. Transverse spectra 

The transverse spectrum of the emerging pions is conveniently expressed as the invariant 
differential yield, Ed 3 N/d 3 p = d 3 'N j 'cPkdy . Due to the boost invariance of the scenario, 
the resulting inclusive yield does not depend on the rapidity y. Moreover, because of the 
rotational invariance in isospace, it is the same for all three pion states. Finally, the overall 
axial symmetry of the geometry guarantees that the resulting inclusive observables have 
azimuthal invariance. Therefore, the transverse spectra depend only on the transverse kinetic 
energy, E k = m k -m 7V . 

In order to have a useful reference for judging the calculated spectra, we introduce the 
equivalent temperature Tbe which is determined by matching the calculated spectrum with 
a corresponding Bose-Einstein equilibrium form. The energy per nucleon determines the 
temperature, 

e = § = %^ , h = [e mfe/TBE - I]" 1 , (19) 

^ V Z^k Jk 

and the appropriate normalization is determined subsequently by the actual number of 
particles N. This fit is made by considering only a certain energy interval, namely 
200 < E k < 1000 MeV. The lower cutoff is made in order to permit the extraction of 
the enhancement of soft pions that is expected to occur. The exclusion of the high momenta 
is made for numerical convenience: the high-momentum modes (which have only a very 
small average occupancy) are not treated with great accuracy in the present study as they 
are computer-demanding but have a negligible effect on the results of present interest. (It 
has of course been checked that an upwards extension of the upper cutoff on Ek, with a 
corresponding better treatment of the high momentum modes, has no import on the results 
extracted.) 

As an example of the outcome of this procedure, the top panel of Fig. shows the 
transverse spectrum extracted from an ensemble of simulation events having To = 240 MeV, 
together with the corresponding Bose-Einstein fit. The fit is seen to be quite good for the 
region E k > 200 MeV within which it has been determined. However, below that region 
the dynamical result is significantly larger than the corresponding extrapolation of the fit. 
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In order to better bring out this enhancement of the soft pions, we prefer to divide the 
calculated spectral profile by the corresponding Bose-Einstein fit, thus getting curves that 
are approximately unity above 200 MeV. The bottom panel of Fig. |5] shows the resulting 
relative transverse spectra for a number of initial temperatures T . It is evident that the 
spectra are to a good approximation of equilibrium form above 200 MeV. (The slight shortfall 
at the higher end of the interval is due to the numerical cutoff in momentum space and it 
has no import on our present discussion, as discussed above.) Furthermore, equally evident, 
there is a significant enhancement of the soft pions over a range of T values, with the effect 
being largest in the transition temperature region where it reaches up to a factor of two for 
the softest pions. (It should be noted that the relatively sudden onset of the enhancement 
below 200 MeV precludes the possibility of obtaining fits of comparable quality by employing 
an adjustable chemical potential.) 

Such a spectral enhancement would probably be within experimental reach. However, 
since the above results were obtained for idealized scenarios of Bjorken matter, they may 
well overestimate the effect. It is therefore important to ascertain to what degree this 
enhancement will persist as an observable signal after the complications of a finite geometry 
have been included. 



III. BJORKEN ROD 

We now turn to the study of a more refined scenario in which the system has a finite 
extension in the transverse plane. Specifically, preserving the longitudinal scaling expansion, 
we shall replace the matter scenario by a rod-like geometry for which the local environment 
changes from that of hot Bjorken matter in the bulk to that of vacuum outside through a 
smooth surface region with a circular cross section. For convenience, we shall refer to such 
a system Bjorken rod. It is illustrated in Fig. |6|. 

The Bjorken rod combines two features exhibited by the systems produced in real nuclear 
collisions: a rapid longitudinal expansion and a finite transverse extension. This latter 
feature implies the presence of a surface region through which the order parameter changes 
from its initially small value in the hot interior to its vacuum value outside. Because it 
is both relatively realistic and reasonably tractable, this model scenario has been widely 
employed in high-energy heavy-ion studies. In the context of DCC studies, it was employed 



already early on by Asakawa et al. [21]. In that work, the chiral field was taken to be 
constant in the longitudinal direction, thus rendering the system effectively two-dimensional 
(with a corresponding reduction of the pressure in the interior). Furthermore, the present 
study employs a more refined initialization procedure which endows the system with more 
realistic features, such as a finite correlation length and a surface region through which the 
field changes smoothly from the hot interior to the exterior vacuum. Finally, whereas the 



work reported in Ref. |2T| was mainly exploratory and rather limited in scope, the present 
study involves quantitative analyzes of the final state in order to make contact with specific 
experimental observables. It should be emphasized that the general conclusions drawn in 
that exploratory work are not affected by the results of the present study. 
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A. Preparation of the rod 



For the rod scenario, we again employ a rectangular lattice in (x, y, rj) and impose periodic 
boundary conditions in all three directions. In the transverse xy plane, the size of the box 
should be sufficient to contain the expanding system at times up to the point when the 
analysis of the final state is made; when this is ensured, the transverse boundary conditions 
are immaterial. In the longitudinal 77 direction, the periodicity approximates conditions 
that would prevail if the cylinder were infinitely long, as in the idealized scenario originally 
considered by Bjorken J34|. Typically the box sides employed are L x = L y = 44 fm, which 



suffices for the systems considered, and L v = 8, corresponding to that many units of rapidity. 

In preparing the system, we wish the local conditions to correspond approximately to 
thermal equilibrium at a temperature that decreases steadily from To inside the bulk region 
to zero in the vacuum outside. The local temperature T has a Saxon- Woods profile, 

T(p) = T [l+e^ ^]- 1 , (20) 

where p 2 = x 2 + y 2 . For the surface width parameter we have used w = 0.8 fm and we show 
results obtained for ensembles of rods with Rq equal to 6 or 10 fm. We expect such values 
of Rq to provide a lower bound on the transverse extensions that may be reached for central 
collisions of gold nuclei at RHIC energies by the time the system has cooled to the specified 
value of Tq. This feature is convenient, since a system with a larger transverse extension 
would be closer to the Bjorken matter scenario considered above (which is approached in 
the limit R — > oo). 

In order to prepare the initial field configuration for the rod, we proceed at first in the 
same manner as for the preparation of the matter scenario addressed above and sample the 
field configuration (4>b ox , V'box) from a thermal ensemble describing macroscopically uniform 
matter within the overall box containing the calculational lattice. This field configuration 
can be uniquely decomposed into its spatial average, the order parameter, and the remainder, 
the fluctuating part of the field, 

0box(P^) = ± + H{P,v) , (21) 
^box(P^) = t + W&V) • (22) 

In order to obtain the field configuration describing the initial state of the rod, (</> rod , V'rod)) 
we rescale the two parts based on the specified local temperature T(p) and then recombine 
them into the desired initial conditions, 

0rod(p^> r o) = #O)[^-0 gs ] +</> gs + Hp) ° [ s 4>(p,v) - 0) +<!>, (23) 

Vw(p,^o) = 9{p)[±-^J +-0g S + Hp)°[W{p,v)-M+±, (24) 

where </> = (/„., 0) and xj) = (0, 0) are the vacuum values. The scaling coefficient for the 
order parameter is obtained from the magnitudes of the corresponding order parameters, 
g(p) = [0q(T(/?)) — / w ]/[0o(^o) — /tt]) where <t>o(T) denotes the magnitude of the order 
parameter in thermal equilibrium at the specified temperature T. The scaling coefficient for 
the fluctuations depends on the particular chiral component, so h is a diagonal 0(4) tensor 
with the following elements, 
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^ - {§y • ^ - (Hi ■ < 25 » 

and o denotes the corresponding 0(4) scalar product. Here ba\ denotes the thermal equi- 
librium value of the fluctuation in the o component of the chiral field and Sit^ is the thermal 
fluctuation in (any) one of the pion components. These are calculated by the semi-classical 
mean-field method |[25|| . (It should perhaps be noted that in principle one needs to take 
account of the fact that generally the order parameter <fi is not fully aligned with the 0(4) 
a direction and hence the a and 7r field fluctuations are not independent. Rather, it is the 
fluctuations parallel and perpendicular to the order parameter, 5<fiu and S(f) ± , that decouple 
(i.e. their quasiparticle mass tensor is diagonal). However, for the quite large systems con- 
sidered here, cf) is practically fully aligned in the a direction and so we do not need to be 
concerned with that distinction in our discussion. Nevertheless, for the sake of generality, 
the numerical treatment performs the proper 0(4) diagonalization and rotation when the 
field preparation is done l2~5fl .) 



By proceeding in the manner described above, we ensure that the local environment, 
as characterized by the order parameter and the field fluctuations, reflects approximately 
thermal equilibrium in matter held at the local temperature T, which decreases from its 
bulk value T to zero according to the prescribed radial profile (|20|) . 

The initial system is illustrated in Fig. [7] for T = 240 MeV. It shows the radial depen- 
dence of a number of key quantities: the local temperature T, the order parameter <p , the 
field fluctuations 5a and 5n, and the effective masses \i a and p^. Since the local temperature 
drops steadily as a function of the transverse distance p, moving out along the abscissa corre- 
sponds to reducing the temperature (though not at a steady rate). As one moves out through 
the surface, the order parameter increases steadily from its reduced value (~ 27 MeV) in the 
hot bulk region towards its vacuum value f n (= 92 MeV), and the local thermal fluctuations 
drop correspondingly towards zero. (Since generally > the fluctuations along a pion 
direction exceed those in the a direction.) The resulting profiles of the effective masses also 
reflect their temperature dependences: Starting from nearly degenerate values (~ 300 MeV) 
in the hot interior, where chiral symmetry is approximately restored, fi n and /i CT diverge 
steadily towards their free values of 138 MeV and 600 MeV, respectively. For higher values 
of the central temperature To, the central value of the order parameter is smaller and the 
effective masses are larger (and even closer in value) and fi a will in fact exhibit a dip in the 
surface region as the local temperature passes through the critical region |25|J36[ . 



B. Rod dynamics 



After the preparation of the field configuration of the rod, the dynamical propagation is 
readily obtained by solving the field equation ([!]) as was done for the matter configurations. 
Since the system is no longer macroscopically uniform, it is more complicated to discuss, 
but it is especially instructive to see how the bulk region develops. For this purpose, we 
probe a circular region in the interior of the rod and extract the corresponding values of the 
order parameter and field fluctuations. The resulting evolution of the order parameter and 
the field fluctuations is illustrated in Fig. [8] for a rod with T = 240 MeV and Rq = 6 fm. 
The corresponding quantities for Bjorken matter prepared with the same value of T are 
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also shown. It is seen that the environment in the interior of the rod evolves in a manner 
quantitatively very similar to that of the corresponding matter scenario throughout the first 
complete oscillation of the the order parameter, which takes about 4 fm/c. At around that 
point in time, the decompression wave has arrived from the rod surface and the self-generated 
transverse expansion now extends over the entire cross section of the rod. As a consequence, 
the relaxation progresses faster than it would in matter where this effect is absent .[] 

This feature is also brought by the the corresponding dynamical paths on the phase 
diagram, as shown in Fig. |]. The initial location of the point for the rod interior differs 
slightly from that of the equilibrium matter because of the finite size of the hollow cylindrical 
volume sampled. But apart from this minor difference, the rod path tracks the matter path 
closely for the first several time units, as brought out by the proximity of the corresponding 
time markers at r = 1,2,3,4,5 fm/c. Beyond that time, the paths begin to diverge, as 
shown by the subsequent time markers at r = 10, 20, 30 fm/c. From around 20 fm/c, or so, 
the rod interior is quite close to the vacuum point, whereas the environment in the matter 
scenario remains significantly excited for a much longer time. 

The quicker relaxation of the field for the rod configuration is also reflected in the behavior 
of the effective pion mass. This is illustrated in Fig. [10] for the same case. Again we see 
that through the first several time units the effective mass in the interior of the rod follows 
closely the evolution of the effective mass extracted for the corresponding matter scenario, 
but it then drops much faster towards the free value. Furthermore, the oscillations in the 
effective mass persist for quite a long time, thus making it possible to achieve a significant 
degree of parametric amplification. 



C. Transverse spectra 

We now consider a number of specific observables that are practically accessible in the 
analysis of actual experimental data. First we consider the transverse spectra of the emerging 
pions, Ed 3 N/d 3 p = d 3 N/d 2 hdy. As was the case for the Bjorken matter scenarios, the 
transverse spectral shape is well approximated by an equilibrium form above kinetic energies 
of 200 MeV or so. We therefore proceed in the same manner and extract a corresponding 
equivalent temperature, Tbe, for each case considered. 

The values obtained for Tbe are shown in Table [I] for a range of initial bulk temperatures 
Tq and a rod radius of R = 6 fm (the values of T B e are rather insensitive to R ) . We note 
that they increase steadily with T , as one would expect. It might at first seem surprising 
that Tbe, the effective temperature characterizing the transverse spectrum of the final pions, 
can exceed T , the specified bulk temperature of the initial state. However, it should be kept 
in mind that the quasipions in the bulk of the initial system have an effective mass fi n that 
exceeds the free value by hundreds of MeV, due to the self-interaction of the thermally 
agitated field. As the system expands and cools, the interaction energy energy is converted 
primarily into kinetic energy and hence the spectrum hardens. 



* The quicker relaxation of the rod configuration is a computational advantage, as it compensates 
somewhat for the larger lattice needed to enable the rod to expand transversally. 
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The resulting ratio between the dynamical spectrum and the corresponding Bose-Einstein 
form is displayed in Fig. 11 for T = 240 MeV, together with the corresponding matter result. 
While the quality of the fit for rods is as good as for matter, the ensuing soft enhancements 
are generally smaller. This is not surprising, since the bulk conditions prevail in only part 
of the rod configuration and, moreover, the transverse expansion tends to suppress the 
oscillatory relaxation causing the enhancement. 

Figure [12] shows the relative excess of soft pions as a function of the initial bulk tem- 
perature T for Bjorken rods with Rq equal to 6 and 10 fm. We see that the effect begins 
to manifest itself at temperatures above 200 MeV, for which the bulk value of the order 
parameter is significantly reduced from its vacuum value, and it exhibits a broad maximum 
at roughly To = 240 — 280 MeV before gradually subsiding at higher To. The larger the 
rod radius Ro, the larger the enhancement, as one would expect, since the system is then 
better approximated by bulk matter. The excess of pions below 200 MeV is of the order 
of ten per cent. Though significantly smaller than that obtained for infinite matter, such 
enhancements are nevertheless still significant and should be well within observable reach. 

Since the excess pions reflected in the spectral enhancement are produced by the coherent 
oscillation of the order parameter, their properties are expected to differ from the regular 
pions. Therefore, in order to better bring this out, it is useful to consider the soft and hard 
pions separately in the following analyses. 



D. Multiplicity distributions 

We now turn to the discussion of the resulting multiplicity distributions. By invoking 



the relation Qlq ), it is straightforward to calculate the expected number of pions in a given 
rapidity interval and within a specified transverse energy bin. Of course, this result varies 
from event to event, since the field configurations vary at the microscopic level as a result 
of the thermal fluctuations in the ensemble of initial states. 

For a given value of the expected number n, the actual number of particles emerging is 
a stochastic variable n, since the state described by a given field has no well-defined particle 
number. Generally, knowledge of the classical field alone does not allow a construction of 
the underlying quantum many-particle state. However, in the present case where we have 
made the common assumption that the final quantum state is of standard coherent form (see 
Eq. (PI1|)), it is elementary to show that the associated multiplicity distribution is of Poisson 
form. Thus, the probability for obtaining a specified number of particles with transverse 
momentum k within a given rapidity interval, n^(yi,y2) = n, is given by 

f) n 

Pn(n) = -r e~ n , (26) 
n\ 

where n = (yi,y2) is the given mean value (|18|). We can therefore readily pick the actual 
number of j-type pions emerging from each particular rapidity slice with the transverse 
momentum k. The Poisson distributions form a closed algebra under convolution, with 
the corresponding mean values being additive, P a * Pb = P a +b- Consequently, the total 
number of particles emerging with a given transverse wave vector k from a given rapidity 
lump can be obtained equivalently by either sampling the multiplicities for each of the 
constituent slices from the corresponding individual Poisson distributions, or by sampling 
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the the total multiplicity from the corresponding combined Poisson distribution. Different 
values of k can be combined in a similar manner, for example to obtain the number of pions 
emerging within a specific energy bin. It is thus clear that the resulting calculated statistical 
multiplicity distributions do not depend on the specific computational procedure employed 
for the stochastic sampling. 



As an illustration, we show in Fig. [L3| both the expected and the actual multiplicity as 
obtained for a single event. It follows from the above remarks that the fluctuations in the 
calculated rapidity density has two different sources. One is the variation in the details of 
the resulting field configuration from event to event and from one rapidity lump to another 



(giving rise to the smooth curves in Fig. [L3|). The other is the subsequent Poisson sampling 
of the actual multiplicity based on the smooth average given by the calculated final field 
(leading to the fluctuating histogram in the figure). 

In order to investigate whether the resulting multiplicity distributions contain deviations 
from the relatively trivial Poisson statistics, we have extracted the associated factorial mo- 
ments. Because of the enhancement of the soft pions seen in the transverse spectra, we 
consider the pions below and above 200 MeV separately. For multi-particle final states, one 
may define the following factorial moments, 

M m = -<N(N-l)---(N-m + l) y . (27) 

Here iV is the number of pions emitted by a given rapidity lump of unit length and the 
average is over all lumps in a sample of events prepared in the same manner (i.e. same bulk 
temperature T and rod radius Rq)\ the contribution vanishes for a lump with a multiplicity 
smaller than the particular order m. For a Poisson multiplicity distribution characterized by 
the mean multiplicity N, the factorial moments are simply given by Ai m = N m . Therfore 
we prefer to consider the reduced factorial moments, A4 m /N m , which are then all unity for 
a Poisson multiplicity distribution. Since N is not available experimentally, we use instead 
the average of the observed multiplicities for each lump, -< iV y. This should provide a good 
approximation, since the lumps are all macroscopically equivalent. 



Figure [14] shows the relative factorial moments obtained for a sample of Bjorken rods 
prepared with To = 250 MeV and Ro = 6 fm. (We consider from now on systems prepared 
with To = 250 MeV, since the soft enhancement is largest near this temperature.) While 
the hard pions appear to be perfectly consistent with pure Poisson statistics, the soft pions 
exhibit a significant non-poissonian behavior. The character of the deviation of the soft 
factorial moments suggests that the source occasionally emits anomalously many pions, as 
one would expect if some modes are especially amplified. This feature is consistent with the 
enhancement of the soft spectrum and it can also be explored by other means of analyzing 
the data for anomalous fluctuations. 



E. Azimuthal emission patterns 

When averaged over events, the emission probability has overall cylindrical symmetry, 
due to the symmetry of the initial ensemble of fields. Similarly, due to the boost invariance of 
the ensemble, the averaged emission probability is independent of the rapidity. And, thirdly, 
the isospin invariance of the ensemble guarantees that all three types of pion behave similarly 
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after the event average has been carried out. Though each individual lump retains these 
symmetries in a macroscopic sense, texture exists at the microscopic level. This important 
feature is manifested as fluctuations in the underlying mean multiplicities n^\yi,y 2 ), from 
one rapidity interval to the next and from event to event. 

The fluctuations in the azimuthal pattern of the asymptotic pion field can be elucidated 
by means of the following multipole moments,^ 

«& } (yi,y 2 ) = ^lE4 j) (yi,y 2 )^e iM ^|- (28) 

The weighting by the transverse kinetic energy Ej. = m^ — m^ ensures that particles emitted 
with very low transverse momentum have little effect, as is reasonable since their particular 
azimuthal direction is physically insignificant. Because of the spectral enhancement of the 
soft pions discussed above, we extract the multipole strength distribution for soft and hard 
pions separately. Finally, E± = ^k-^fc denotes the total mean transverse kinetic energy of 
the soft or hard pions emitted by a given lump. The division by E± ensures that ck is always 
unity. (It should perhaps be added that the rather large transverse dimensions of the spatial 
lattice ensures that the corresponding resolution in the wave number k is sufficiently fine 
to render the extracted multipole moments coefficients numerically reliable for values of the 
multipolarity M comfortably beyond the point when their values have become insignificant.) 

In each individual event, and for each of the three pion types separately, the multipole 
moments are calculated for the pions arising from each separate rapidity lump. Because 
of the independence of different events, the isospin symmetry, and the longitudinal boost 
invariance, respectively, all the individual lumps present equivalent sources and hence we 
may perform a subsequent average, yielding -< a M y. 

The resulting ensemble average azimuthal multipole strength function is displayed in 
Fig. [15] for the soft modes of the asymptotic pion field. Two different Bjorken rod scenarios 
are displayed, both having the initial bulk temperature is To = 250 MeV, but the initial 
radius Rq being either 6 fm or 10 fm. In order to get an idea of the significance of the ex- 
tracted values, we have calculated an auxiliary set of moments, <S^ ndom , obtained by turning 
each individual emission direction </>k by a random angle before evaluating the multipole 
moment. 

The azimuthal coefficients for soft pions, considered as a function of the multipolarity 
M, start out well above the noise level, then exhibit a gradual drop-off. For high multipo- 
larities, the moments approach zero because the high-frequency components of the field are 
suppressed by the initial thermal weight. 

The behavior of the multipole strength provides some information on the domain struc- 
ture of the source, as we shall now discuss. First we note that the attenuation of the strength 
with increasing multipolarity M can be reasonably well approximated by a simple gaussian 



^The second-order multipole coefficient gives a rough indication of the quadrupole moment of the 
emission pattern and is related to the coefficient v% introduced in Ref. [37] for the discussion of 
elliptic flow. 
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fit.0 The fit makes it possible to extract the dispersion of the curve, AM. Employing the 
usual inverse relationship between conjugate variables, we may then obtain the azimuthal 
correlation length, A0 = 1/ AM. Thus, if the effective transverse size of the source is R, 
the spatial domain size in the source is given by K-RAcf). For the two scenarios having the 
initial rod radius Rq equal to 6 and 10 fm, we extract AM as 5.03 and 8.84, leading to values 
of A(p equal to 11?4 and 6?5, respectively. Since RoA<ft is thus approximately the same in 
the two scenarios, the spatial correlation length at the effective emission time appears to be 
independent of the initial radius of the rod. This is as one would expect if that quantity 
reflects the conditions in the source at the time of effective decoupling. 

The corresponding multipole coefficients extracted for the hard pions are shown in 
Fig. [16|. They are initially an order of magnitude smaller than those for the soft pions 
but they exhibit a significantly gentler fall-off with the multipolarity M, suggesting a cor- 
respondingly narrower angular correlation. This is consistent with the general expectation 
that harder pions exhibit a smaller spatial correlation length. As for the soft modes, the 
wider rod yields smaller multipole coefficients. 

The above discussion has been based on the pion field, rather than the resulting particles. 
The field is a relatively smoothly varying function of position and it governs the mean 
number of quanta in the various elementary modes. The actual final state of particles is 
then determined stochastically, as described in the preceding. This feature tends to add a 
certain degree of noise to the azimuthal pattern exhibited by the underlying field, as we shall 
now illustrate. When the analysis is made in terms of the emitted particles, as is necessarily 
the case in actual experiments, the definition of the multipole moments should be modified 
appropriately, 

*m = ^-|£^e w *»| , (29) 

where n enumerates the N particles emerging in the particular kinematic domain under 
consideration and E n = m n — m n is the transverse kinetic energy of an individual particle 
emitted in the azimuthal direction <p n . Furthermore, E± = Y^n^ni is the total transverse 
energy of the N particles, so that a is unity. 

Figure IT7| shows the multipole strength coefficients o.m extracted for the case with Rq = 



6 fm. As was the case for the coefficients extracted from the underlying field, the particle- 
based coefficients acquire significant values for low values of the multipolarity M. As the 
multipolarity is increased, aM approaches the statistical noise level, a^| ndom , obtained as 
above by randomizing the direction of each emitted pion. This is what one would expect 
since the angular distribution is now a sum of ^-functions and thus contains all the high 
multipolarities equally. To better bring out the features of the non-statistical behavior, we 
subtract the average value of the noise, (a random ), before making the fit (hence some of the 
resulting values are slightly negative for large M). 



Hn the numerical treatment, the transverse pion momenta are given on a cartesian grid which has 
a four-fold C4 symmetry with respect to azimuthal rotations. Consequently, those moments whose 
multipolarity M is divisible by four exhibit anomalous strength, as is best seen for the high values 
of M . We have therefore omitted those multipolarities from the fitting procedure. 



14 



Once the average noise has been subtracted, the multipole strength distribution based 
on the emitted particles is quantitatively very similar to the one based on the underlying 
asymptotic field. This is true for both values of Rq shown. It thus appears possible to gain 
experimental access to the correlation properties of the field itself, even though the particles, 
on which the observations are necessarily based, are obtained from the field by a stochastic 
process. 

The above analysis suggest that an extraction of azimuthal moments for observed soft 
pions may yield non-trivial information on the domain structure in the system. In this 
undertaking, it would probably be useful to explore the effect of adjusting the cuts in both 
rapidity and transverse energy, so as to tune the extraction procedure to the specific character 
of the sources produced. 



F. Neutral pion fraction 

The distribution of the neutral pion fraction / = N^o/N^ has received considerable at- 
tention as a DCC diagnostic because an idealized fully isospin-polarized source would yield 
an anomalously wide distribution, Po(f) = l/2y/J, rather than a narrow peak around the 
average value / = |. Indeed, experiments have been devoted to the search for anomalous 
behavior of the neutral pion fraction |T3|JT^] . However, such an undertaking is quite diffi- 



cult, partly because the signal is expected to be carried by only the soft pions and partly 
because the extraction of / requires the simultaneous event-by-event measurement of both 
charged and neutral pions. These practical problems not withstanding, P(f) remains a very 
instructive quantity in the DCC context and we have therefore extracted it as well. 

As above, we consider soft and hard pions separately. In each case, we calculate the 
neutral pion fraction / for each rapidity lump and construct its distribution P(f) by binning 
the / values obtained from all lumps and events. f\ Since all extracted distributions P(f) 
appear well-peaked around their common average of one third, it suffices for our discussion 
to characterize a given distribution by its variance er| = (f 2 ) — (f) 2 . It is elementary to 
show that the extreme idealized distribution Po(f) has the variance a 2 = 4/45 and that 
N such sources of equal strength have the variance a 2 = <Jq/N P5[| . More generally, one 
may consider an assembly of N different sources, with individual distributions Pi(fi). If the 
strength (i.e. the total multiplicity) of a given source is Ui, the neutral fraction for the entire 
assembly is given by / = Yh u ifi/ Si v % an d its distribution is therefore given by 

- JV - "I N N 

P(f) = UJdfiWi) «(/-EWE"i)- ( 3 °) 

_i=l J i=l i=l 

It is easy to verify that the mean value of / remains equal to one third. Moreover, a 
slightly more elaborate calculation yields the expression for the variance of the combined 



§ Since isospin symmetry guarantees that that the labeling of the three cartesian isospin axes may 
be permuted without any physical effect, the sampling error on P(f) can be reduced by making 
the corresponding three entries into the histogram from each individual source, rather than merely 
a single entry [29|. 
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distribution, 



a 



2 S« ,J ; _2 



(31) 



If the isospin polarizations of all the elementary modes k are totally independent, then the 
above formula fl3"T|) applies and it is possible to predict the resulting width 07 by substituting 
the mean multiplicities = J2j ^-k f° r the strengths z/j. However, in actuality the pion field 
changes smoothly from one mode to the next and so the isospin directions of neighboring 
modes are somewhat correlated. Consequently, we expect the actual 07 to be significantly 
larger than the lower bound provided by (1311) and indeed this turns out to be the case. 

In the present discussion, we calculate the neutral fraction both on the basis of the 
asymptotic fields, which govern the mean multiplicities for each elementary mode, and on the 
basis of the resulting particles, which are picked stochastically from those mean multiplicities 
(see above). Table [Ij lists the resulting widths 07 for Bjorken rods with radii Rq equal to 
either 6 or 10 fm. We notice that when the transverse size of the system is increased, the 
distribution becomes narrower, as expected since more domains can be accommodated across 
the source. In addition to considering the usual lumps covering one unit of rapidity each, we 
have also calculated 07 by separating each lump into its five constituent slices. Naturally, 
the associated decrease of the source size causes P(f) to broaden. However, the the fact that 
the variance grows by less than a factor of five is a reflection of the fact that the effective 
domain size extends beyond a single slice, so that neighboring slices are not independent 
with regard to their isospin alignment. Finally, we notice that the distributions based on 
the particles (which are the observable quantities) are generally significantly broader than 
those based on the underlying fields. 

It is also instructive to study how P(f) evolves when the data is separated into azimuthal 
sections. Thus, each individual lump is subdivided into m equally large azimuthal sections 



and P(f) is calculated for each such source. The result is illustrated in Fig. [TSj which shows 
the width of P(f) as a function of the number of angular sections m, for both soft and hard 
pions. Generally P(f) broadens as the width of the angular section is decreased, as one 
would expect since the effective number of contributing domains goes down. However, 07 
will not approach the value for a single idealized source, as m — > 00, because the transverse 
extension of the source exceeds the correlation length and so the source will never become 
fully aligned in isospace. This feature is similar to what was obtained when idealized rods 
were sliced ever finer in rapidity |38|. Moreover, as we have seen above, the soft pions lead 



to a significantly broader distribution than the hard ones. 



IV. DISCUSSION 

In the present work, we have employed the standard linear a model to study systems 
prepared with a cylindrical geometry and endowed with a longitudinal Bjorken scaling ex- 
pansion. Such a scenario is an idealized representation of what might be produced in a 
central ultrarelativistic nuclear collision. After their preparation with an interior tempera- 
ture in or above the chiral transition region, the system is left free to evolve self-consistent ly, 
until the combined longitudinal and transverse expansion has brough it into the simple 
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asymptotic regime of free dispersion where a well-defined analysis can be made. We have 
studied a number of observables related to the emerging pions, for the purpose of exploring 
the persistence of signals reflecting the preceding non-equilibrium relaxation of the chiral 
field. 

Most important, and also most easily checked experimentally^ is the emergence of a 
surplus of pions with transverse kinetic energies below 200 MeV, relative to the expectation 
based on a simple extrapolation of the harder part of the spectral profile, which turns out to 
be quite well approximated by a Bose-Einstein form. The excess of soft pions calculated in 
this manner depends on the initial bulk temperature of the rod, T , as well as on its radius, 
Rq, and it is largest for T = 240 — 280 MeV. In the optimal range of T , the total relative 
excess of pions below 200 MeV amounts to about 10% and 15% for R equal to 6 and 10 
fm, respectively. Since the transverse radius of the source is likely to be larger than 10 fm 
for central collisions at RHIC, the present rod calculations are expected to underestimate 
the effect. Thus, if indeed the effect is present, it ought to be readily visible in the data. 

Since the additional soft pions arise from the oscillatory relaxation of the chiral order 
parameter in the interior of the source, they are expected to be non-statistical and partially 
isospin polarized. In order to elucidate the degree to which remnants of these features are 
identifiable in the final state, we have performed various event-by-event analyses on the 
calculated results based on samples of several tens of separate rod evolutions. 

In one type of analysis, we have extracted the factorial moments of the pion multiplicity 
distribution associated with a rapidity interval of unit width. We found that while the hard 
pions display a perfect Poisson multiplicity distribution, the soft pions show a significant 
non-poissonian behavior, with the reduced factorial moments growing approximately linearly 
with order. A similar multiplicity analysis should be readily possible with the multi-particle 
data collected by the STAR TPC, for example. 

In another series of analyses, we sought to elucidate the azimuthal structure of the 
final emission pattern by extracting the corresponding multipole strength function. These 
analyses were made both for the asymptotic field, which is not directly detectable, and for 
the resulting particles. The strength was found to decrease steadily with multipolarity in an 
approximately gaussian manner, thus making it possible to extract an angular correlation 
length. These studies suggested that the spatial correlation length in the source (the "domain 
size") is independent of the initial rod radius, but significantly larger for the soft pions. 
Moreover, these features were not affected by the stochastic transition from field to particles, 
thus suggesting that such angular analysis of the actual data may be informative. 

Finally, we extracted the distribution of the neutral pion fraction for the event samples 
considered, since this quantity has played a central role in the discussion of disoriented chiral 
condensates (though it is practically difficult to measure). Information about the correlation 
length, and hence the domain structure, was obtained by dividing each rapidity lump into 
azimuthal segments before calculating the neutral fraction. We found that the fraction 
distribution widens steadily as the number of segments is increased and, importantly, that 
the width associated with the soft pions is several times larger than that for the hard pions. 



** It should be noted that the spectral profile is a one-particle observable and thus easier to 
address experimentally than quantities requiring true event-by-event analysis. 
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Even though the present study is more refined than earlier treatments, with regard to the 
dynamical scenario and the analysis of the outcome, it is still highly idealized in a number of 
respects. Especially important are the following limitations. First, at the most basic level, 
the linear a model is an effective field theory containing merely a few mesonic fields, rather 
than the basic quark-gluon degrees of freedom. This simplification may become especially 
important at the early stage when the high excitation is expected to largely dissolve the 
hadronic correlations in the interior of the system. Second, even at the effective level, a 
full non-equilibrium real-time quantum-field treatment is presently far beyond practicality. 
Instead, we are resorting to a semi-classical treatment. While perhaps quite adequate in 
many respects |4(| , this treatment is expected to generally underestimate the enhancement in 
the pion yield caused by the non-equilibrium dynamics |35| . Third, the commonly employed 
SU(2) version of the linear a model entirely ignores strangeness. However, this degree of 
freedom is expected to become significantly agitated once the temperature exceeds the mass 
of the strange quark, as is the case at the early stage of the evolution. The extension 
of the model to SU(3) has been found to have significant effects [|39[| . Fourth, the initial 



state is obtained by sampling a (suitably modulated) field configuration from a thermal 
ensemble. The sampling procedure relies on the mean-field approximation and thus the 
initial state contains only a minimum of many-particle correlations. This feature might 
become reflected in a corresponding paucity of non-trivial correlations in the final state. 
An analogous problem concerns the transcription of the final classical field configuration 
to a specific quantum state. Since our prescription assigns a standard coherent state, it 
is likely to provide an oversimplification of the many-body correlations present. (Indeed, 
recent studies suggest that the parametric amplification mechanism behind the spectral 
enhancement in fact destroys the simple coherent form of the state |^5|.) Fifth, the present 
study addresses only part of the system. A more realistic scenario would contain a halo 
of hadrons outside the hot cylinder and their interaction with the emerging pions might 
well be significant (although the fact that the interesting signals are carried by slow pions 
may give the surrounding material time to disperse before having any significant effect). 
In addition, of course, the cylindrical geometry is merely an idealization of actual collision 
geometries occurring. (In addition to these theoretical simplifications , our study ignores 
the deficiencies inherent in the detection system.) 

Because of these various approximations and idealizations, our calculations should not 
be considered as directly comparable with experimental data. Rather, our study serves to 
elucidate the importance of some of the complications that are inevitably present in any 
real system produced by high-energy collisions, especially the rapid longitudinal expansion 
and the presence of a surface region between the agitated interior and the tranquil vacuum 
outside. Generally speaking, our results suggest that while one should not expect spectacular 
DCC signals in the actual data, carefully designed quantitative analysis might bring out any 
anomalous DCC-type behavior in the data. In particular, the observables considered in our 
analyses of the simulation results, especially the transverse spectrum, the factorial moments, 
and the multipole strength, may be useful in this regard. 
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TABLES 



TABLE I. The fitted BE temperatures. 



T (MeV) 


200 


225 


250 


275 


300 


T BE (MeV) 


256 


284 


324 


367 


412 



The table shows the equivalent temperatures Tbe obtained by fitting a Bose-Einstein profile to 
the calculated transverse pion spectra in the kinetic-energy interval between 200 and 1000 MeV, 
for Bjorken rods with an initial bulk temperature of To and a radius of R® = 6 fm. Rods with 
i?o = 10 fm lead to very similar values. 



TABLE II. Width of the neutral pion fraction distribution. 



Source type 


Fields 


Particles 


Source size 


Lumps 


Slices 


Lumps 


Slices 


Soft (Rq=6) 


5.34 


6.51 


11.74 


26.98 


Hard (Rq=6) 


1.29 


1.82 


8.15 


19.54 


Soft (#0=10) 


3.25 


4.14 


6.66 


14.16 


Hard (R =10) 


0.80 


1.08 


4.70 


10.65 



The table shows the dispersion aj of the neutral pion fraction distribution P(f) (given in per cent), 
obtained for two ensembles of Bjorken rods having an initial bulk temperature of To = 250 MeV 
and a transverse radius Rq equal to either 6 or 10 fm. The neutral fraction has been calculated 
using either the mean multiplicities (given by the field strength) or the emerging particles (obtained 
stochastically based on the given mean). The sources are taken to be either the lumps covering 
one unit of rapidity or their constituent slices which each cover Ay = 0.2. Pions with transverse 
kinetic energies below (soft) and above (hard) 200 MeV have been considered separately. 



20 



FIGURES 




FIG. 1. Time evolution of the field. 
The time dependence of tfio, the magnitude of the order parameter 0o (top panel), and 5(f), the 
associated dispersion of the field fluctuations (bottom panel), in Bjorken matter (a longitudinally 
expanding box) prepared at various initial temperatures, Tq, as indicated. 
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FIG. 2. Phase evolution. 
The time evolution of the expanding box is projected onto the chiral phase diagram, in which 
the abscissa is the order parameter 0o (the magnitude of the average value of the chiral field within 
the box) and the ordinate is the filed fluctuation 5(j) (the dispersion in the field around its average 
value), for the four different initial temperatures To. The equilibrium path is shown by the upper 
dotted curve while the lower dotted curve delineates the region of instability within which the field 
is supercritical. 
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FIG. 3. Effective mass in Bjorken matter. 
The time evolution of the square of the effective pion mass, /J%(t) (divided by the square of its 
free mass m n ) for Bjorken matter prepared at four different initial temperatures To at the proper 
time tq = 1 fm/c. 
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FIG. 4. Enhancement factor. 



The number enhancement coefficient Xjq obtained as described in Ref. |35| on the basis of the 
time dependent effective masses given in Fig. ||. The coefficient Xn expresses the relative increase 
in the population of a given pion mode, as characterized by the momentum p = Tike. 
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FIG. 5. Transverse spectrum for Bjorken matter. 
The top panel shows the final transverse pion spectrum, d 3 N/d 2 pj_dy, and the associated 
equilibrium spectrum obtained by fitting the result for the longitudinally expanding Bjorken box 
with a Bose-Einstein form within the energy interval 200-1000 MeV, for the initial temperature 
To = 240 MeV. The bottom panel shows the ratio between the dynamically generated spectrum 
and the fitted equilibrium form, for four different initial temperatures Tq. 



25 




FIG. 6. The initial geometry of the Bjorken rod. 
The initial field configuration describes a rod-like system subject to a longitudinal scaling 
expansion of the Bjorken type, for which the local boost rapidity is given by y = zjr. 
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Initial transverse profile of Bjorken rod 




Transverse distance p (fm) 



FIG. 7. Initial rod profiles. 
The transverse profiles of the local temperature T, the order parameter <po, the mean fluctuation 
of the field in a given 0(4) direction, 5a and Sir, and the corresponding effective masses ii a and 
/iyr, for a rod with a radius of Rq = 6 fm and a central temperature of Tq = 240 MeV. 
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FIG. 8. Time evolution of the field. 
The time dependence of (f>o, the magnitude of the order parameter (increasing curves), and 
Scj), the associated dispersion of the field fluctuations (decreasing curves), in both Bjorken matter 
prepared with To = 240 MeV (solid curves) and the corresponding Bjorken rod with Rq = 6 fm 
(dashed curves). The information for the rod has been obtained by averaging over a hollow cylin- 
drical volume, 1 fm < p < 3 fm. (The region very close to the symmetry axis has been excluded in 
order to avoid the enhanced fluctuations in the center caused by the strict axial symmetry of the 
geometry.) 
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FIG. 9. Phase evolution of the rod. 
The phase evolution of the interior of the Bjorken rod prepared with a bulk temperature 
of To = 240 MeV and a radius of Ro = 6 fm (dashed path, filled diamonds), compared with 
the corresponding evolution of Bjorken matter shown in Fig. ^ (solid path, open circles). The 
equilibrium path is shown by the short-dashed curve while the dotted curve delineates the region 
of instability within which the field is supercritical. The information for the rod has been obtained 
by averaging over a hollow cylindrical volume, 1 fm < p < 3 fm. The symbols give the location of 
the paths at successive proper times r = 1, 2, 3, 4, 5, 10, 20, 30 fm/c. 
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FIG. 10. Effective mass. 
The time evolution of the square of the effective pion mass, /J%(t) (divided by the square of its 
free mass m n ) for a Bjorken rod (dashed curve) prepared with To = 240 MeV and Rq = 6 fm as well 
as the corresponding result for Bjorken matter taken from Fig. [| (solid curve). The information 
for the rod has been obtained by averaging over a hollow cylindrical volume, 1 fm < p < 3 fm. 
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FIG. 11. Spectral shape for Bjorken rods. 
The ratio between the final transverse pion spectrum, d 3 N/d 2 pj_dy, and the associated equi- 
librium spectrum obtained by fitting the dynamical result with a Bose-Einstein form within the 
energy interval 200-1000 MeV, for a longitudinally expanding rod (dashed curve) having an ini- 
tial radius of Ro = 6 fm and with an initial bulk temperature of To = 240 MeV, as well as the 
corresponding result for Bjorken matter (from Fig. |5|). 
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FIG. 12. Excess of soft pions for Bjorken rods. 
The relative excess of pions for Bjorken rods with initial radii Rq equal to either 6 or 10 fm, 
as a function of the initial bulk temperature To. The excess has been obtained by subtracting the 
yield of pions with kinetic energies below 200 MeV from the corresponding equilibrium spectrum 
obtained by fitting the dynamical result with a Bose-Einstein form within 200-1000 MeV. 
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FIG. 13. Rapidity density. 
The rapidity density of pions for a single event prepared with To = 250 MeV, as obtained by 
dividing the number in each slice by its width Ay = 0.2. The expected multiplicity based on the 
calculated asymptotic field is indicated by the smooth solid curve, while the histogram shows the 
actual number as picked from the corresponding Poisson distributions. The top curves include all 
three pion types, whereas the bottom curves show the neutral pions only. 
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FIG. 14. Reduced factorial moments. 
The reduced factorial moments M m /N rn obtained for a sample of rods prepared with Rq = 6 fm 
and To = 250 MeV displayed as a function of the order to, for either soft (dots) or hard (squares) 
pions. The error bars arise from the sampling over all lumps and events and the points have been 
slightly offset for clarity. 



34 



Azimuthal pattern of the soft modes 



o.o 



T =250 MeV, R =6 fm 



■ Multipole strength 
Gaussian fit 
Randomization 




Bjorken rod 
T =250 MeV, R =10fm 

• Multipole strength 
- Gaussian fit 
Randomization 



o^o ooooo " 



12 16 20 

Azimuthal multipolarity M 



24 



28 



FIG. 15. Multipole strength for soft modes. 
The azimuthal multipole strength distribution (in per cent) for the soft modes of the asymptotic 
pion field, a (solid circles), extracted for an ensemble of events starting with an initial rod radius 
Rq of either 6 or 10 fm and an initial bulk temperature of To = 250 MeV. The solid curve shows 
the corresponding gaussian fit. The open diamonds indicate the noise level d^ ndom , obtained by 
randomizing the azimuthal direction of each pion mode before evaluating the multipole moment. 



35 



0.20 



Azimuthal pattern of the hard modes 
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FIG. 16. Multipole strength for hard modes. 
The azimuthal multipole strength (in per cent) of the hard modes of the asymptotic pion field, 
a, extracted for an ensemble of events starting with an initial rod radius i?o of either 6 or 10 fm and 
an initial bulk temperature of To = 250 MeV, together with the corresponding gaussian fit. The 
open diamonds indicate the noise level d^ ndom , obtained by randomizing the azimuthal direction 
of each pion mode before evaluating the multipole moment. 
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Azimuthal pattern of the soft particles 
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FIG. 17. Multipole strength for soft particles. 
The azimuthal moments (in per cent) extracted for an ensemble of events starting with an 
initial rod radius Ro = 10 fm and an initial bulk temperature of To = 250 MeV. These moments 
have been extracted from the many-particle final states obtained by sampling the actual particle 
numbers from the underlying asymptotic fields used in Fig. 15. The noise level averaged over the 
multipolarity M, (a random ), has been subtracted from the displayed values (they amount to 2.48 
and 0.81 per cent for Rq equal to 6 and 10 fm, respectively) and the solid curve shows the gaussian 
fit to the resulting difference. 
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FIG. 18. Neutral pion fraction. 
The dispersion crj of the neutral pion fraction distribution P(f) for sources that have been 
obtained by dividing each rapidity lump into m equally large azimuthal sections of width 2n/m. 
An ensemble of rods with a initial bulk temperature To = 250 MeV and radius Ro = 10 fm. The 
results for both soft and hard pions are shown. These moments have been extracted from the 
many-particle final states resulting from the underlying asymptotic fields used in Fig. [l5|. 
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